Minimal two-sphere model of the generation of fluid flow at low Reynolds numbers 
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Locomotion and generation of flow at low Reynolds number are subject to severe limitations due 
to the irrelevance of inertia: the "scallop theorem" requires that the system have at least two degrees 
of freedom, which move in non-reciprocal fashion, i.e. breaking time-reversal symmetry. We show 
here that a minimal model consisting of just two spheres driven by harmonic potentials is capable of 
generating flow. In this pump system the two degrees of freedom are the mean and relative positions 
of the two spheres. We have performed and compared analytical predictions, numerical simulation 
and experiments, showing that a a time-reversible drive is sufficient to induce flow. 
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Microscopic systems capable of generating flow are 
very common in Nature, and may prove inspirational for 
bio- mimetic micro- and nano-pumps and swimmers 
Assemblies of motile cilia are found in various eukaryotic 
living systems. In humans, for example, they cover the 
epithelial tissue of important organs, including the lungs, 
ventricles in the brain, and the oviduct in the female 
reproductive apparatus 0]. They transport fluid along 
their surface in a given direction by controlling the ef- 
fective drag coefficient to change between a "power" and 
a "recovery" stroke Q. This is one simple way to sat- 
isfy the "scallop theorem" Q that sets very strict phys- 
ical requirements for swimming and pumping at small 
velocity to viscosity ratio at the microscale, where the 
Reynolds number is to good approximation zero. A nec- 
essary condition, in order to pump, is that the sequence 
of the system's configurations has to break time-reversal 
symmetry The scallop theorem applies to pumps as 
well as swimmers [f|, so no net flow will occur unless 
the generating motion is non-reciprocal. This implies a 
minimum of two degrees of freedom, with which time- 
reversal symmetry can be broken by an appropriate se- 
quence of moves M, 0] • The design of micro and nano- 
fluidic devices 0, Q which mimic biological examples is 
an emergent field of research with potential applications 
in medicine and biotechnology [10| . If future nano-bot 
swimmers and pumps might be made through a process 
of self assembly, the question is how few components are 
necessary to generate flow, and how simple can the sys- 
tem be. 

In this paper we describe a minimal model of pump, 
inspired by the three-sphere swimmer [f| @, fill [l2j where 
the actuated motion along one axis reduces the descrip- 
tion to a one dimension [13[ . The two-sphere system re- 
duces further its complexity. In the limit of low driving 
frequency and for average bead separation larger than 
the bead diameter, the hydrodynamic interaction is de- 
scribed by the Oseen tensor 0, [HI, an d the equations 
of motion are simple enough to allow for explicit calcu- 
lations. The analytical results back up numerical calcu- 
lation and experimental data to confirm the surprising 
result that two beads subject to harmonic potentials can 



generate a net flow even when the external drive is recip- 
rocal. 




FIG. 1. (Color online) The two-bead micropump setup. Bead 
L is held in a stationary optical trap, while bead R is sub- 
ject to a time-reversible driving force from a switching opti- 
cal trap. This assembly is capable of generating flow. The 
harmonic trap potential is shown schematically overlaid on 
snapshots taken from the experiments. 



The pump is composed of two spheres of radius a la- 
belled with L (Left) and R (Right), positioned on the 
x-axis at an initial distance d as illustrated in FIG. [T] 
Each sphere is subject to an harmonic potential, which 
is realized experimentally by an optical trap, anchored 
to the laboratory reference frame. Bead R is actively 
driven: its confining potential switches between two po- 
sitions along the x direction separated by a distance e, 
with a frequency l/2r. The position of the minimum 
varies as a function of time as a a square wave accord- 
ing to a(t) = x[2nr, (2n + l)r](f), where x denotes the 
indicator function and n is a positive integer. Bead L is 
in a stationary potential. There are two distinct phases, 
corresponding to the values a = 1 , of the active drive, 
which constitute a basic cycle of dynamics. 

Despite the fact that the trap movement is reciprocal in 
time, the external actuation by means of springs cannot 



2 



guarantee an instantaneous balance of the active forces, 
so the pump is not instantaneously force- free [l|, [l5[ . Un- 
like for swimmers, however, this is not problematic: a 
pump is a spatially confined system and its center of 
mass lies within a bounded region. In this system this 
circumstance is also a necessary condition, assuring that 
two degrees of freedom are accessible for pumping. No- 
tably, a weaker condition holds: the pump is force-free 
on average over a cycle period. 

In a low Reynolds number fluid, the actively driven 
particle undergoes a purely dissipative dynamics, just like 
a forced oscillator that relaxes towards the minimum of 
its confining potential. While bead R moves, it inter- 
acts hydrodynamically with bead L, which can vary its 
position thanks to the softness of the harmonic poten- 
tial. Their dynamics is described, in the regime where 
a <C d, using the Oseen tensor approximation (l6j |. De- 
noting with xl the coordinate of L and similarly for R, 
the governing equations read 



XL 



(1) 



xr = — cj( xr, — (d — ecr) ) — uj— xl', 



where we have introduced the three parameters 7 = 
6iTr]a, u> = fc/7, A = 3a/2 and the relative distance 
r = xr — xl- Here rj is the fluid viscosity, 7 is the 
Stokes' drag coefficient and k is the stiffness of the spring. 
The equations show that the geometric parameters of the 
model have a clear interpretation: the inverse of d sets 
the strength of the hydrodynamic coupling between the 
spheres, and e is the oscillation amplitude of the actively 
driven particle. 

The model reveals an intriguing property. The pres- 
ence of two active drives, without any constraint, would 
provide the system with two "obvious" degrees of free- 
dom. However, the temporal dependence of the active 
mechanism described here is symmetric under time re- 
versal, as can be seen by inspecting FIG. Q] from top to 
bottom and vice- versa. Thus, at first sight it might ap- 
pear that the system cannot generate net flow. Instead, 
the left-right symmetry is broken. 

The pumping can be quantified by focusing on the bead 
in the resting trap, L; an asymmetry in its trajectory 
reveals an asymmetry in the flow field. We define the 
order parameter 
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to quantify the magnitude and the direction of the flow 
as a function of the parameters d, e and r. Its physi- 
cal interpretation is as follows. Imagine that bead L is 
an isolated sphere attached to a spring, and subject to 
the same flow field as the one generated by the pump. 
Then Hooke's law gives F = kAx, allowing to measure 
the equivalent mean force exerted by the pump. Using 



Stokes' law this flow field can be converted into a mean 
velocity of the fluid v — F/7. 

The analysis of the two phases of motion helps to un- 
derstand how the symmetry breaking occurs [? ]. I) in 
Phase I, the hydrodynamic coupling between beads L 
and R has a strength of the order of 1 jd and increases as 
the spheres approach to reach their minimum distance. 
According to the positions of bead L, the fluid is pushed 
in the sense x < 0. Eventually bead L is restored back to 
equilibrium xl = 0, and in this relaxation some fluid is 
dragged in the opposite direction. II) in Phase2, hydro- 
dynamic coupling is stronger [on the order of l/(d — e)]. 
The dynamics looks similar to the mirror image of the 
previous one, with bead R dragging bead L in the sense 
of x > 0, but now the stronger coupling moves a greater 
amount of fluid. The overall effect is a net thrust in 
this direction. An example of such motion for bead L 
is illustrated in FIG. EJa). The mismatch in hydrody- 
namic coupling between end of Phase 1 and beginning 
of Phase 2 is the root of the symmetry breaking and is 
made possible by the softness of the driving potentials. 
Such a phenomenon is analogous to the soft swimming 
described in [I?} . As we see the pumping direction is de- 
termined by the position of the actuated particle: if bead 
R is active, the pumping is in the direction of x > and 
vice versa if bead L is active. 

Introducing the reduced distance u = r — d and mean 
coordinate c = (xl+xr) /2, in the approximation of small 
oscillations the equations of motion can be expanded as 
power series in the parameter u/d. With this change, 
equations ((T|) decouple into an equation for u, indepen- 
dent of c, and a linear equation for c involving both vari- 
ables. The dynamics at the leading order in u is obtained 
by substituting l/(it + d) ~ 1/d in the resulting equa- 
tions. This corresponds to the study of the linearized 
system and a careful analysis shows that no pumping is 
achieved. Physically, this fact can be understood by in- 
terpreting 1/r as an effective drag felt by the center of 
mass c. When r is approximatively constant, then the 
drag doesn't change during the two phases of dynamics 
causing therefore no net thrust on c and thus on the fluid. 

Pumping arises as a non-linear effect which can be seen 
already at the next order of expansion. According to 
l/(d + u) w 1/d — u/d 2 , the reduced distance u has to 
satisfy of a set of Riccati equations depending on the 
parameter a: 



u = P 
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where we have defined V a := —ecrui (1 — \/d), Q a := 
-oj (1 - \/d + a{\e)/d 2 ) and K := -uiX/d 2 . The center 
of mass equation maintains its linear character, 
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for the reduced variable c a := c— (d — ea) /2. Both equa- 
tions can be solved for a = 0, 1. Further one must impose 
appropriate conditions on the solutions, representing: (i) 
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the continuity of the whole solution in the middle of the 
cycle where u (t) = iti(r) and Cq(t) = Ci(t); and (ii) 
the steady state condition for which the positions at the 
beginning of the cycle coincide with those at the end, 
given by uo(2t) = iti(O) and co(2t) = ci(0). Finally, 
using the inverse transformation from u,c to get xl and 
xr and taking the temporal average of xl, we find that 
the order parameter Ax shows a net pumping over two 
phases of dynamics. 

FIG. [2] reports the plots of the resulting expression as 
functions of e, t, d for typical experimental values of a 
and u. There, the geometrical variables e and d follow the 
straightforward trend where larger oscillation amplitudes 
e generate larger flow. Also, larger distance d is related 
to smaller generated flow. However, the plots show also a 
non-trivial behavior as a function of the temporal param- 
eter r. The form of the solution suggests a natural way 
to rescale the parameters, defining the nondimensional 
quantities £* := e/a, d* := d/a and t* := tlo. There are 
two different regimes, i) for small values of r* the flow 
increases and ii) for large values of r* it decreases. The 
reason for the latter is that in this regime the beads have 
sufficient time to relax in the harmonic traps, so that the 
area described by the function xl reaches a maximum 
value, after which it remains constant at increasing r*. 
Dividing this area by r* gives a decreasing function of r* . 
Interestingly, at the intersection of these two regimes the 
plot shows that Ax presents a maximum, thus defining 
an optimal pumping region. Due to the particular form 
of the function Ax, however, this has to be evaluated nu- 
merically. For the typical experimental values in use it 
corresponds to r w 165ms. 

For t> 1/w, where we can compare with our experi- 
mental data, the expression of Ax simplifies considerably. 
At the leading order the result is a simple power law de- 
pendence 

a ~ 16t* (d*) 3 ' K ' 

Integration of the full equations of motion ([1]) is also 
performed numerically by means of Taylor's method [llj]. 
Comparison with the analytical results shows a good 
agreement, despite the low-order expansion of the analyt- 
ical solution, as can be seen by looking at FIG. GHa),(b). 
This is not surprising as the Oseen tensor is a descrip- 
tion valid for of large separation of the beads, and in this 
regime the perturbative result is a good approximation 
to the exact solution. 

The dynamics has been investigated experimentally 
by means of an optical tweezer, described fully in [7j. 
A laser beam (IPG Photonics, PYL-1-1064-LP, A = 
1064nm, Pmax = 1.1W) is focused through a water im- 
mersion objective (Zeiss, Achroplan IR 63x/0.90 W), 
trapping from below. The laser beam is steered via 
a pair of acousto-optic deflectors (AA Opto-Electronic, 
AA.DTS.XY-250@1064nm) controlled by custom built 
electronics, allowing multiple trap generation by time 
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FIG. 2. (Color online) Analytical characterization of the 
pump, (a) Aa; as function of e and r. The trend as a func- 
tion of £ is intuitive: the larger the oscillation amplitude, the 
larger the generated flow. There is a non trivial behavior 
as function of r, showing that the maximum pumping is ob- 
tained at maximum e for some value of r in the middle of the 
scale, d here is fixed to the experimental value of 6/im. (b) 
Aa; as function of d and r. We see the same phenomenology 
as in (a). Here the effect of the distance is to decrease the 
pumping monotonically. (c) Plot of Aa; as a function of e and 
d. Again, it shows a monotonic dependence from the distance 
and e. Looking at these plots, is easy to understand that the 
scaling exponents of Aa; are functions of all these parameters. 
In particular, the scaling law of Aa; with e depends on r. In 
all these cases the values of uj and a are set to typical values 
taken from the experiments. 



sharing, with sub-nanometer position resolution. Instru- 
ment control and data acquisition (70 frames per second, 
with an exposure time of 13 ms) are performed by custom 
software. 

A typical experiment consists of trapping two silica 
beads (3.0/im diameter, Bangs Labs) in a solution of glyc- 
erol (Fisher, Analysis Grade) and is divided in two parts: 
in the first calibration stage all the traps are kept at rest, 
and the beads undergo only Brownian motion confined 
by the traps. During the second stage we reproduce the 
cycle of FIG.JT] with lasers traps, iterating the sequence 
many times [7j. The whole procedure lasts typically 6 
minutes, during which we collect about 40 000 frames. 
We performed 3 runs for each set of parameters. 

By analyzing images by correlation filtering and two 
dimensional fitting, we obtain the beads' position with 
subpixel (around lnm) resolution. The expected values 
of Ax are below the limit of the experimental resolution 
(indeed the smallest forces measured here correspond to 
O.lnm displacements), and therefore we characterize the 
pumping indirectly, relying on the measurable asymme- 
try of peaks of bead L's position 0- Regarding each cycle 
as an independent realization of a stochastic process, we 
construct the mean dynamic cycle of L, see FIG. [3Ja). 
Compared to the equilibrium position of the bead, that 
can be determined from the mean cycle itself but only for 
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FIG. 3. (Color online) Experiment, analytical calculations 
and simulations agree in quantifying flow generation by the 
two-bead micropump. (a) As bead R is driven by the optical 
trap, it causes displacements of bead L around its equilibrium 
position. Markers show the position of bead L averaged over 
many cycles. The analytical solution, (green dashed curve) 
and simulation data (blue continue curve) closely reproduce 
the experiments. The displacement values p+ and p_ at the 
peaks are an accurately measurable quantity. Matching the 
simulations to the peak values enables the observed p+ and 
P- to be converted to force by a one-to-one mapping; (b) The 
mean force exerted by the pump on the fluid is measured from 
the motion of the bead in the stationary trap, and it grows as 
a function of the pump stroke length e. Markers and lines are 
as in panel (a). In addition the approximate solution, Eq. (JS} 
is shown (black dash dotted line). 

large values of r* when the dynamics is fully relaxed, we 
indicate the maximum with p + and the minimum with 
p_ . We define the algebraic sum S := p+ +p~ to quantify 
the asymmetry of motion and convert it to Aa; with the 
aid of simulations. This procedure allows to compare ex- 
tremely small forces, of the order of 5 x 10~ 4 pN, which, 
to our knowledge, are the smallest forces measured with 
optical traps. 

In FIG.[3£b) we plot the corresponding mean force ob- 
tained from the experiments at varying e and compare 



this result with analytical predictions and simulations. 
The experimental values considered are a = 1.5/xm, 
d = 6/xm, t = 640 ms, w = 0.022±0.001(ms)- 1 , T = 25° 
C and a trap stiffness value of k = 5.32 ± 0.71 pN//xm. 
The results show a good agreement between the mea- 
sured force and the predicted values. Approximate ana- 
lytical solution ([3]) gives a good description of pumping. 

It is interesting to compare the effectiveness of the cur- 
rent minimal pump to the related three-bead model sys- 
tem Q. Close to the optimal pumping region, and for 
matching values of the stroke e and inter-particle distance 
d, the average force exerted on the fluid by the three- 
bead model exceeds the two-bead companion by about 
one order of magnitude. The poorer performance with 
two beads is not surprising, and minimality is obtained 
at the expense of performance. However we would like 
to point out that the two pumps have also a profoundly 
different nature. In the three-bead model, pumping is 
achieved by moving the lateral beads in a non-reciprocal 
fashion. The direction of the flow is determined by the 
first moving trap and can be reverted inverting the trap 
moves. The nature of the drive in the two-bead model 
prevents all this and the pumping direction is uniquely 
determined by the disposition of the active trap, as dis- 
cussed above. 

In conclusion, an extremely simple system composed 
of just two spherical beads, only one of which is actuated 
by a time-reversible trap movement, is shown to be ca- 
pable of generating flow. A key property of the system 
is that the beads are held and driven by soft potentials. 
This allows the two-bead system to explore two degrees 
of freedom, thus satisfying the "scallop theorem" . The 
simplicity of this elementary pump makes it possible to 
understand the fluid dynamics analytically, and suggests 
this as a feasible micro-pump that could be deployed ex- 
perimentally in the context of microfluidic systems. 
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